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Non-perturbative renormalisation with domain wall fermions 
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We present results from a study of the renormalisation of both quark bilinear and four-quark operators for 
the domain wall fermion action, using the non-perturbative renormalisation technique of the Rome-Southampton 
group. These results are from a quenched simulation, on a 16^ X 32 lattice, with /3 = 6.0 and Ls = 16. 



1. Introduction 

Domain wall fermions are increasingly be- 
ing used for many studies of physically interesting 
quantities. It is therefore important to have reli- 
able values for the renormalisation coefficients of 
the operators needed in these calculations. Here 
we will discuss the application of the nonpertur- 
bative renormalisation technique of the Rome- 
Southampton group [Q , the Regularisation Inde- 
pendant (RI) scheme, to quark bilinears and four- 
quark operators using the domain wall fermion 
action. Rather than concentrating on the numer- 
ical values for the renormalisation coefficients, the 
emphasis will be on highlighting tests of the chiral 
properties of domain wall fermions and on inter- 
esting aspects of the non-perturbative renormali- 
sation technique. 

All the results presented are from a study using 
a lattice with dimensions 16^ x 32, an extent in 
the fifth dimension of 16 points, a domain wall 
height (M5) of 1.8 and a Wilson gauge action with 
/3 = 6.0. 

2. The RI scheme 

A full description of the RI scheme may be 
found in jl], but the relevant details will be 
quickly summarised here. 

The matrix element of the operator of interest. 
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between external quark states, at high momenta, 
in a fixed gauge (in this case Landau gauge was 
used) is calculated. The external legs of this ma- 
trix element are then amputated. A renormalisa- 
tion condition may be defined by requiring that, 
for the renormalised operator at given scale, a 
chosen spin and colour projection of this quan- 
tity is equal to its free case value. 

3. Quark BiUnears 

For the case of flavour non-singlet fermion bi- 
linear operators, we write 



Zr [uTid] 



(1) 



with Ti = {l,ltj.,l5,lfil5,o-fj.u}, where i repre- 
sents whatever indices the gamma matrices have. 
Taking projections of the amputated matrix ele- 
ments with the same gamma matrices as in the 
definition of the operator, normalised such that 
this trace is unity in the free theory, gives the 
renormalisation condition 



with 



Ar 



(2) 



niron=0 



Tr [r,r, 



-TT[r,{[uT,d]uip)dip)}AMp] ■ (3) 



Invariance of the action under axial flavour 
transformations requires that Za = Zy and Z5 = 
Zp. The extent to which these relations hold for 
domain wall fermions is therefore an excellent test 
of their chiral properties. Figure |l| shows data for 
A^ — Ay versus (ap)^ for several values of the 
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input mass, mf. While there is a small differ- 
ence between Aa and Ay at low momenta, this 
difference becomes smaller as either the momen- 
tum becomes larger or the mass becomes smaller, 
suggesting that the bulk of this difference is due 
to the effects of chiral symmetry being sponta- 
neously broken by the vacuum and softly bro- 
ken by the mass respectively. Within the statis- 
tics there is therefore no signal for a difference 
between Za and Zv- Figure shows A,9, Ap 
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Figure 1. A plot of Aa — Ay versus (ap)^, showing 
that there is no significant difference between Za 
and Zv- 

and their ratio after the mass poles have 
been subtracted by fitting the mass dependence. 
Again, within the accuracy of the method there 
is no signal for a splitting of Zs and Zp. 

4. Four Quark Operators 

The renormalisation of four-quark operators is 
much more complicated than that of the quark bi- 
linears considered above. Mixing both among the 
four-quark operators and between these operators 
and lower dimensional operators must be consid- 
ered. However, these operators are of much in- 
terest phenomenologically, in particular because 
of their relevance for the calculation of Bk and 
matrix elements of the AS* — 1 Hamiltonian [Q. 

4.1. A5 = 2 Hamiltonian 

The AS* = 2 Hamiltonian, which is relevant for 
calculating Bk, is one of the simplest four-quark 



Figure 2. A plot of As, Ap and their ratio, af- 
ter the mass poles have been subtracted, show- 
ing that the any splitting between Zs and Zp is 
smaller that the resolution of this study. 

renormalisation problems, as mixing with lower 
dimensional operators is not possible, ft consists 
of the operators Oab = sTAdsTB d, where Ta 
and r B are arbitrary gamma matrices except for 
the condition that there be no free indices. 

As with the bilinears, renormalisation condi- 
tions may be placed on these operators by consid- 
ering their amputated matrix elements between 
external quark states. As mixing between four- 
quark operators is possible, a single renormali- 
sation condition no longer suffices; however, by 
varying the projectors used, all the renormalisa- 
tion coefficients can be fixed . 

To completely specify the renormalisation con- 
dition, the momentum configuration for the ex- 
ternal quark states must also be chosen. For 
such calculations to be useful, however, we are 
required to use the same momenta configuration 
used in existing perturbative matching calcula- 
tions, which take all the momenta to be equal. 

The renormalisation structure of these opera- 
tors is strongly constrained by chiral symmetry 
Previously presented results |^ have con- 
firmed that the parity conserving part of the oper- 
ator (1 — 75) dsj^ (1 — 75) d does not signifi- 
cantly mix with any of the four wrong chirality 
operators it could possibly mix with. Chiral sym- 
metry also predicts that the renormalisation fac- 
tor for the parity conserving and parity breaking 
components of the operator are equal. The ra- 
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tio of these factors is shown in Figure H, with no 
evidence for chiral symmetry breaking visible for 
moderately high momenta. 




Figure 3. The ratio of the parity positive and 
parity negative renormalisation constants for the 
AS = 2 Hamiltonian. No evidence of explicit 
chiral symmetry breaking can be seen. 

Another prediction of chiral symmetry is that 
the operators = sdsdi 575^5756? do not mix 
with one another. The element of the inverse mix- 
ing matrix resulting from a naive application of 
the RI scheme renormalisation conditions to this 
case is shown in Figure land, as can be seen, this 




Figure 4. The element of the mixing matrix con- 
necting O'^ and , which chiral symmetry con- 
strains to be zero. 

is significantly non-zero. The mass behaviour is 



also very distinctive, with the signal being fur- 
ther away from zero the smaller m/ is. This is 
not, however, a signal for a large mixing between 
these two operators, but an interesting systematic 
effect in the matrix elements we are calculating. 

A simple understanding of this effect can be 
gained by noting that while we are interested in 
the matrix elements of operators between quark 
states, we are still evaluating this matrix element 
in a theory with propagating pseudo-goldstone 
bosons. We may approximately represent the ef- 
fects of the pseudo-goldstone bosons by an effec- 
tive interaction of the form 



C {n d-y^u + Wf^d tt} , 



(4) 



with each such interaction between the quark 
states and the pseudo-goldstone bosons giving 
rise to a propagator of the form 



1 



{p - qf + Ml 



(5) 



where p and q are the momenta on the incoming 
and outgoing quarks. 

In the case we are considering, the problem is 
that all the external quark states have equal mo- 
menta and any interaction between them through 
the vertex in Equation ^ will give a contribution 
to the matrix element oc = l/mquark- 

Figure |^ shows the same quantity as Figure IJ, 
for nif = 0.02, except that the momenta of the 
incoming external states, p, and outgoing states, 

can be different. The data on the graph is such 
that for each point p^ — q^ with four distinct 
blocks of values of p^ being shown, while (p — qY 
varies between the different points. For all of the 
points except the four that are distinctly lower, 
p ^ q, while for those four points (jp — q)^ = 0. 
The points for which {p — q)'^ ^ are clearly 
suppressed as is expected from Equation ^. 

4.2. AS = 1 Hamiltonian 

The mixing problem for the AS* = 1 Hamil- 
tonian is much more complicated that for the 
AS" — 2 Hamiltonian, both in the number of four- 
quark operators that may mix and because of the 
presence of mixing with lower dimensional oper- 
ators. The latter must be subtracted before the 
standard RI scheme conditions are applied, and 



trapolate exactly to zero at the same point as the 
renormalised mass (m = — TOi-qsI^), using a Hnear 
extrapolation, the constant piece is small and is 
of the order of potential subleading effects. 



Figure 5. The same quantity as shown in Figure ^ 
except that the momenta for the external states 
is such that the pseudo-goldstone boson pole is 
suppressed for all but four points on this graph. 

it is the feasibility of this subtraction for domain 
wall fermions that we will touch on here. 

The fact that the RI conditions are applied 
in a fixed gauge and for off-shell quark states, 
makes the problem of mixing with lower dimen- 
sional operators very difficult |Q , with many op- 
erators that must be subtracted. However, this 
number is greatly reduced if chiral symmetry is 
respected. Also, when calculating the renormal- 
isation coefficients in the context of a practical 
study, perturbation theory must be relied upon 
for the running of the effective Hamiltonian and 
scheme matching. Operators that do not appear 
in these calculations (because they are higher or- 
der in perturbation theory) need not be consid- 
ered unless they are power divergent in the lattice 
spacing. 

As it is chiral symmetry that makes perform- 
ing this subtraction feasible, here we will give 
one example of a subtraction coefficient and its 
chiral properties. The AS = 1 operator Og = 
{sadp)L J2qi'ip1o')R may mix with the lower di- 
mensional operator sd. If chiral symmetry is 
not respected, then the needed subtraction co- 
efficient is cubically power divergent in the lat- 
tice spacing with only subleading mass depen- 
dence. Chiral symmetry, however, constrains 
the leading behaviour of this coefficient to be 
oc (I/a^)(ms-|-md). Figure || shows the behaviour 
of this subtraction coefficient versus m f , for a sin- 
gle momentum value. While this does not ex- 
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Figure 6. The coefficient of sd for its subtraction 
from Oq. The strong mass dependence is a sign 
that, to a good approximation, chiral symmetry 
is respected. 

5. Conclusions 

In conclusion, for the lattice parameters 
used in this study, domain wall fermions have 
excellent chiral properties. Also, the non- 
perturbative renormalisation technique of the 
Rome-Southampton group seems to work well as 
long as one is careful when choosing the momen- 
tum configuration used. 

REFERENCES 

1. G. Martinelh et al, NPB 445, 81 (1995). 

2. Y. Shamir, NPB 409, 90 (1993); Y. Shamir 
and V. Furm an, NPB 439, 54 (1995). 

3. D. Becirevic |hep-lat/9909082| ; CudeU et al, 
Phys. Lett. B454 (1999), 105. 

4. T. Blum (RBC), R. Mawhinney (RBC), J. 
Noaki (CP-PACS) these proceedings. 

5. A. Donini et at, Eur.Phys.J CIO, 121 (1999). 

6. C. Dawson (RBC), Nucl. Phys. Proc. Suppl. 
82, 854-856 (2000). 

7. C. Dawson et a l, NPB 514, 313 (1998). 

8. T. Blum et al, |hcp-lat/0007038| . 



